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LONG-TERM  GOALS 

The  long-term  goals  are  to  understand  surf  zone  processes  related  to  rip  current  systems  through  field 
observations.  Rip  currents  occur  commonly  on  most  beaches  and  dominate  many.  It  is  recognized  that 
beaches  with  straight  and  parallel  contours  are  not  a  stable  morphologic  configuration  whereas  more 
complex  beaches,  which  support  the  existence  of  rip  current  morphology,  are  stable  and  more 
common. 

OBJECTIVES 

The  research  objectives  of  the  proposed  work  focus  on  obtaining  new  observations  of  the  three- 
dimensional  structure  of  the  rip  current  system  utilizing  a  suite  of  in  situ  instruments.  In  addition,  a 
fleet  of  30  inexpensive  surfzone  drifters  were  constructed  and  deployed  to  evaluate  Fagrangian  rip 
current  associated  observations.  The  second  related  effort  applies  a  numerical  model  (Delft3D)  to 
evaluate  the  dynamics  of  the  rip  current  system  and  its  interaction  with  the  surface  wave  field  and 
bottom  topography.  These  new  observations  will  be  used  to  validate  Delft3D  and  extend  our 
understanding  of  rip  current  processes. 

APPROACH 

We  (MacMahan,  Stanton,  Reniers,  Thornton,  Gallagher,  Brown,  Brown,  Henriquez,  Stockel,  Cowen, 
Wycoff,  and  Morrison)  conducted  a  Rip  Current  Experiment,  RCEX,  at  Sand  City,  Monterey  Bay,  CA 
in  April-May  2007.  A  combination  of  in  situ  Eulerian  measurements,  remote  sensing  techniques,  and 
Lagrangian  measurements  were  deployed.  The  Eulerian  measurements  consisted  of  two  primary 
arrays:  1)  a  cross-shore  array  of  co-located  pressure  and  digital  electromagnetic  current  meters  (PUV) 
and  ADCPs  along  the  axis  of  the  rip  channel,  and  2)  an  alongshore  array  of  PUVs  and  Paroscientific 
pressure  sensors  (PARO;  Figure  1).  The  two-axis  current  and  pressure  PUV  sensors  provide  data  for 
analyzing  alongshore  wavenumber- frequency  spectra,  and  a  cross-shore  array  of  bottom-mounted 
ADCP’s  captures  the  cross-shore  variability  in  vertical  flow  structure  within  a  rip  channel.  Remote 
sensing  systems  consisted  of  video  and  HE  radar. 
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30  surfzone  drifters  with  accurate  GPS-tracking  were  deployed  for  three  hours  for  seven  different  days 
under  varying  wave  and  tidal  conditions  to  quantify  the  spatial  variation  in  mean  Lagrangian  flow, 
dispersion,  and  diffusion.  The  GPSs  after  post-processing  have  an  absolute  position  error  of  <  0.3  m 
and  speed  errors  of  <  3  cm/s  (MacMahan  et  ah,  2008).  High  resolution  velocity  measurements  over 
large  area  are  required  to  map  the  complete  cell  circulation  of  rip  currents. 

Concurrent  numerical  model  predictions  of  the  local  hydrodynamic  conditions  were  performed  to  help 
in  the  deployment  of  the  surfzone  drifters  and  the  execution  of  jet-ski  surveys.  These  computations  are 
based  on  the  transformation  of  deep  water  directional  spectra  to  the  nearshore,  including  the  wave 
groups  to  simulate  the  three  dimensional  infragravity  time  scale  surfzone  circulations. 

Jeff  Brown  (University  of  Delaware  graduate  student),  Rob  Wyland  (NPS  tech),  Ron  Cowen  (NPS 
tech),  Jim  Lambert  (NPS  Tech),  Jon  Morrison  (NPS  student),  Jamie  MacMahan  (NPS)  constructed  the 
surfzone  drifters.  MacMahan,  Stanton,  Renters  (RSMAS),  Thornton(NPS),  Gallagher  (Franklin  and 
Marshall),  Brown,  Brown,  Henriquez  (Delft),  Stockel  (NPS  tech),  Cowen,  Wycoff  (NPS  tech),  and 
Morrison  were  responsible  for  instrument  deployment,  maintenance,  data  archiving,  removal,  drifter 
depoyments,  and  bathymetric  surveys.  Currently,  there  are  5  students  utilizing  the  dataset  for  Masters 
thesis.  Two  students  (Jeff  Brown  and  Jenna  Brown)  are  from  the  University  of  Delaware  and  three 
students  (Jon  Morrison,  and  Andrea  O’Neil)  from  the  Naval  Postgraduate  School.  In  addition  Martijn 
Henriquez  (University  of  Delft)  will  use  the  data  obtained  from  the  near  bed  observations  of  velocity 
and  sediment  fluxes  as  part  for  his  PhD  thesis  to  examine  intra-wave  sediment  transport  processes. 

WORK  COMPLETED 

We  successfully  deployed  an  alongshore  array  of  PUVs  and  a  cross-shore  of  ADCPs  within  a  rip 
channel.  We  performed  7  drifter  deployments  under  various  wave,  tidal,  and  surfzone  flow  conditions. 
We  performed  5  bathymetric  surveys  over  the  course  of  the  experiment.  Delft3D  has  been  coupled  to 
the  global  wave  model  Wavewatch  III,  where  the  fireq.  directional  spectra  from  Wavewatch  III  are 
used  as  input  to  the  embedded  SWAN  model  within  Delft3D  to  calculate  the  transformation  of  deep 
water  wave  conditions  to  the  nearshore  zone.  The  SWAN-derived  nearshore  wave  conditions  were 
subsequently  used  as  input  for  the  wave  group  modeling  suite  within  Delft3D  resulting  in  predictions 
of  the  three  dimensional  surfzone  circulation  patterns  on  the  infragravity  time  scales  and  longer.  Next 
the  combined  three  dimensional  Lagrangian/Eulerian  flow  field  was  extracted  from  the  Delft3D  flow 
computations  and  used  to  predict  the  dispersion  of  surfzone  drifters.  The  GPS-drifter  observations 
have  been  analyzed  to  yields  estimates  of  surfzone  diffusion  and  dispersions  using  two-particle  and 
cluster  statistics.  Post-experiment  modeling  has  been  performed  to  verify  modeling  concepts  and 
numerical  implementations  in  Delft3D,  using  the  observed  frequency  directional  spectra  at  12.8  m 
water  depth  to  construct  the  offshore  wave  boundary  (Renters  et  ah,  2008)  and  compare  with  GPS- 
drifter  inferred  surface  velocities  and  Eulerian  in  situ  measured  sub-surface  velocities. 

RESULTS 

The  experiment  was  highly  successful  in  providing  many  necessary  observations  for  evaluating  rip 
current  dynamics. 
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Lagrangian  Observations 


The  rectified  time-averaged  breaking  wave  patterns  for  four  representative  drifter  observation  patterns 
are  plotted  in  Figure  1  (left  column).  Drifter  observations  are  subdivided  into  three  primary  flow 
patterns:  1)  symmetric  and  asymmetric  rip  currents  that  were  observed  when  Hmo.ism  >  1  m,  2) 
wandering  currents  that  was  a  combination  of  rip  current  flow  pattern  and  sinuous  alongshore  flow 
pattern  that  were  observed  when  Hmo,i3m>  Im,  and  3)  sinuous  alongshore  currents  that  were  observed 
when  Hmo,i3m<  Im. 

For  the  first  time,  a  synoptic  circulation  pattern  of  a  rip  current  is  observed  in  the  field  with  enough 
resolution  to  calculate  the  vorticity  field  (Figure  2a, b).  Two  different  rip  current  circulation  patterns 
were  observed:  1)  symmetric  rip  current  circulation  (Figure  2a)  and  2)  asymmetric  rip  current 
circulation  (Figure  2b).  The  symmetric  rip  current  has  onshore  flows  over  the  shore-connected  shoals 
that  transition  to  alongshore  flows  (feeder  currents)  near  the  shoreline  which  converge  and  move 
offshore  in  the  rip  channel,  and  transition  to  alongshore  flows  that  diverge  near  the  surf  zone  edge 
inducing  a  counter  (positive  vorticity)-  and  a  clockwise  (negative  vorticity)  rotational  element  (rip 
current  vortex).  These  two  opposing  rotational  elements  are  a  rip  current  system  that  is  generally 
bounded  in  the  cross-shore  by  the  surfzone  width  and  in  the  alongshore  by  the  "quasi"  alongshore 
periodic  rip  channels.  The  faster  velocities  occur  near  the  outer  edges  of  the  rotational  elements,  with 
the  slowest  velocities  occurring  near  the  vortex  centers.  Asymmetric  rip  currents  are  similar  to 
symmetric  rip  currents,  except  that  the  vortex  occupies  one  rip  channel  and  one  shore-connected  shoal 
and  has  up-  or  down  coast  preference.  Areas  of  opposing  vorticity  occur  near  the  edges  of  adjacent  rip 
current  cells.  Maximum  absolute  vorticity  occurs  in  the  center  of  the  rip  current  cells.  The  rip  current 
pattern  and  surfzone  vorticity  estimates  are  consistent  with  numerical  model  estimates  [Renters 
et  ah,  2007;  amongst  others]. 

Lagrangian  rip  current  circulation  patterns  are  quantitatively  documented  in  the  field  for  the  first  time. 
These  Lagrangian  observations  provide  useful  insights  about  the  behavior  of  rip  currents  and  surfzone 
exchange.  Surprisingly,  the  rip  current  flow  field  consisted  of  semi-enclosed  vortices  that  retained 
material  within  the  vortex  center  and  were  contained  within  the  surf  zone.  The  drifter  observations 
increase  toward  the  center  of  the  rip  current  vortices  for  energetic  conditions  when  the  circulation 
pattern  is  coherently  organized.  Talbot  and  Bate  [1987]  found  surfzone  diatom  blooms  were  correlated 
with  larger  breaking  waves  and  rip  current  events,  and  the  largest  concentrations  of  surfzone  diatoms 
were  toward  the  center  of  rip  current  vortices.  Our  drifter  results  support  their  findings  and  indicate 
that  despite  the  large  velocities  in  the  rip  channel;  diatoms,  pollution,  and  other  fioatsom  are  not 
necessarily  transported  offshore  but  are  instead  retained  and  concentrated  within  the  rip  current 
vortices  inside  of  the  surf  zone.  Unlike  ocean  vortices  that  trap  and  transport  material,  the  rip  current 
vortex  is  fixed  in  one  location,  as  it  is  coupled  to  the  surfzone  morphology,  and  only  retains  the 
material  and  does  not  transport  it. 

The  observed  semi-enclosed  surf  zone  circulation  patterns  have  significant  importance  to  Navy 
personnel  with  regard  to  beach  safety.  For  most  rip  current  drifter  deployments,  swimmers  were  often 
caught  in  the  rip  current  while  releasing  the  drifters.  Swimmers  would  float  and  ride  the  circulation 
pattern  back  to  wading  depths  on  the  shore-connected  shoal,  mimicking  the  drifter  response.  The 
drifters  can  be  considered  proxies  for  floating  humans.  For  the  asymmetric  rip  current  day,  GPSs  were 
placed  on  a  few  swimmers  as  they  freely  floated  in  the  rip  current  system.  An  example  of  one  person 
who  floated  around  a  rip  current  for  two  revolutions  is  shown  in  Figure  2.  The  average  swimmer 
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revolution  was  approximately  7  minutes,  similar  to  average  drifter  revolution  [Brown  et  al.,  2008]  and 
dye  revolutions  [MacMahan  et  al.,  2008].  In  general,  beach-goers  are  told  to  "swim  parallel"  or  along 
the  shoreline  when  pulled  into  a  rip  current  (www.usla.org).  The  alongshore  currents  feeding  the  rip 
current  can  be  as  large  as  the  rip  current  (Figure  2).  These  results  suggest  that  if  the  swimmer  remains 
calm  and  afloat,  they  will  be  returned  to  wading  depths  in  less  than  three  minutes.  However,  there  are 
times,  though  infrequent  10%,  that  they  may  exit  the  surf  zone.  No  swimmer  exited  the  surf  zone 
during  the  experiment. 
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Figure  1.  left  column)  Three  video  cameras  overlooking  the  experiment  site  were  mounted  on  a 
tower  18.2  m  above  the  mean  sea  level,  (middle  column)  Drifter  position  and  speed  (colorline)  tracks 
for  the  four  representative  flow  patterns.  The  colorbar  at  the  top  represents  drifter  speed,  (right 
column)  For  drifter  observations  described  herein,  the  biases  in  sample  density  and  temporal 

intervals  [Davis,  1991] 

Swenson  and  Niiler,  1996]  are  minimized,  as  drifters  repeatedly  looped  around  a  rip  current  rotational 
element  and/or  drifters  were  repeatedly  re-released  providing  temporally-evolving  iterative 
observations  in  space.  Owing  to  these  factors,  each  drifter  was  considered  an  independent  observation 
when  in  a  particular  defined  binned  area.  If  the  drifter  re-entered  the  same  area,  it  was  considered  a 
new  independent  observation,  but  only  for  t>lg/U  has  elapsed,  where  Ig  is  the  length  of  the  bin  and  U 
is  the  average  speed  in  the  box  for  all  drifter  observations.  Five  independent  observations  are 
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considered  as  statistieally  reliable  [Spy dell  et  al.,  2007].  Bathymetric  contours  are  plotted  in  the 
baekground  in  gray.  Colorbar  at  the  top  represents  the  number  drifter  observations  for  10m  by  10m 
bins  with  greater  than  four  independent  observations. 

Surfzone  Dispersion  and  Dijfusivity  on  a  Rip  Channeled  Beach 

Absolute  and  relative  diffusivity  are  eomputed  within  the  surf  zone  using  30  position- traeking  drifters 
released  in  elusters  (4-12  drifters)  deployed  on  7  days  with  different  wave  forcing  and  tidal  elevations 
on  a  sandy  beach  with  quasi-periodic  alongshore  spaeed,  0(125  m),  rip  ehannels  (eurrents)  at  Sand 
City,  Monterey  Bay,  CA.  Diffusivity  and  dispersion  were  found  to  be  larger  on  days  with  rip  eurrent 
flow  patterns  and  larger  waves  and  are  eomparable  to  previous  surf  zone  estimates  at  the  asymptotie 
limit.  The  eross-shore  absolute  diffusivity  asymptotes  measured  here  (0.9  -  2.2  m  /s)  on  days  with 
rip  eurrent  flow  patterns  are  similar  in  magnitude  to  earlier  dye  measurements  on  a  beaeh  with  rip 
currents  (0.08  -  5.9  m  /s),  but  alongshore  absolute  diffusivities  measured  here  (2.8  -  3.9  m  /s)  are  an 
order  of  magnitude  larger  than  the  0.03  -  0.2  m  /s  measured  by  Inman  et  ah,  [1971].  Absolute  diffusion 
in  the  rip  eurrent  flow  patterns  is  initially  dominated  by  the  eross-shore,  as  found  by  Inman  et  al. 
[1971],  but  for  large  t  the  diffusion  becomes  alongshore  dominated  similar  to  Johnson  and  Pattiaratehi 
[2004].  The  large  initial  diffusivity  followed  by  a  sharp  decline  results  from  the  expansion  and 
contraction,  as  the  particles  are  quiekly  dispersed  initially  then  begin  to  re-group,  slowing  the  spread. 
The  rip  current  causes  material  to  diffuse  quiekly  for  t<90  s  in  the  eross-shore  {k^x  =  4.9  -  6.1  m  /s), 
before  deereasing  to  an  asymptotie  oseillation  (Kxx  =  0.9  -  2.2  m  /s),  while  alongshore  material 
diffusion  is  initially  (t<170  s)  smaller  than  eross-shore  diffusion  and  asymptotes  at  a  larger  value  (%;  = 
2.8  -  3.9  m^/s).  The  eross-  and  alongshore  absolute  diffusivity  modulates  at  -300  s  eorresponding  to 
the  average  eireulation  time  for  a  rip  eurrent.  For  the  sinuous  alongshore  flow  field,  whieh  more 
elosely  represents  the  eonditions  observed  by  Spy  dell  et  al.  [2007],  the  alongshore  diffusivities 
measured  here  are  similar  in  magnitude  (2.4  -  2.7  m  /s)  to  theirs  (2.0  -  4.5  m  /s).  Initial  eross-shore 
diffusivities  for  the  alongshore  current  systems  (2.6  -  3.7  m  /s)  are  larger  than  those  observed  by 
Spydell  et  al.  [2007]  (0.7  -  1.5  m  /s)  though  they  deerease  to  more  eomparable  values  at  later  times 
(0.1  -  2.5).  On  days  exhibiting  the  same  flow  pattern,  absolute  diffusivities  are  smaller  for  days  with 
smaller  waves  and  shorter  periods,  resembling  Spydell  et  al.  [2007],  though  the  exaet  relationship 
between  flow  pattern,  wave  eonditions  and  diffusivity  eannot  be  generalized  due  to  the  laek  of  data. 
The  eross-shore  diffusivity  on  yearday  139  exemplifies  the  differenee  bathymetrie  variations  and  flow 
pattern  ean  make  on  the  diffusivity.  Wave  heights  and  period  are  similar  to  day  1  of  Spydell  et  al. 
[2007],  but  the  sinuous  alongshore  flow  restricts  particle  dispersion  to  the  surf  zone  and  the  cross¬ 
shore  diffusivity  falls  to  zero,  while  for  Spydell  et  al.  [2007]  the  eross-shore  diffusivity  reaehes  a  non¬ 
zero  asymptote. 
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Figure  2.  a,b)  Rip  current  velocity  data  are  sorted  into  10m  by  10m  bins  and  averaged  over  the 
deployment  duration  (~3  hours)  for  that  day  for  bins  with  greater  than  four  independent 
observations  (Figure  1,  right  column).  The  red  arrows  in  the  upper  right-hand  corner  provide  vector 
scales.  The  local  bottom  morphology  is  contoured  in  the  background.  Vorticity  is  calculated  in 
discrete  form  using  a  weighted  central  difference  method  and  plotted  behind  the  velocity  vectors 
(left  column  -  a,b).  Note  since  bin  quantity  is  computed  from  nine  bins  surrounding  the  center  bin, 
there  are  a  small  number  of  observations  that  occur  without  a  velocity  observation  within  that 
center  bin.  Owing  to  the  statistical  uncertainity  of  the  spatially-weighted  central-difference 
kinematic  estimate,  a  Monte  Carlo  simulation  is  performed  using  the  spatial  mean  velocity  field  and 
randomizing  the  velocity  standard  deviation  for  each  grid  cell.  The  mean  standard  deviation  for 
Monte  Carlo  simulations  is  0.019s'^,  which  we  consider  as  our  error.  The  hot  and  cold  vorticity 
estimates  are  generally  larger  than  the  error,  c)  The  circles  on  the  dotted  line  represent  minutes 
starting  at  zero  (green)  and  ending  at  14  (red)  for  the  human  drifter  track.  Two  revolutions  of 
human  track  are  plotted  (magenta  for  the  F‘  revolution  and  blue  for  the  2'"^  revolution).  The  local 

bottom  morphology  is  contoured  in  the  background. 


Battjes  [1975],  assuming  straight  and  parallel  contours  used  an  energy  balance  (- 


dEc 


gx 


dx 


=  )  to  solve 


for  the  turbulent  eddy  diffusion  by  assuming  the  production  of  turbulence  generated  is  locally  equal  to 
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dissipation  by  breaking  waves,  Sz,,  where  E  is  wave  energy  and  the  group  veloeity  in  the  x  direetion 
Cgx~  Cg  for  the  assumed  normally  ineident  waves.  The  turbulent  eddy  viseosity  given  by: 


=  aH 


\P) 


(1) 


where  p  is  density  of  seawater,  with  the  seale  of  the  turbulenee  given  by  H,  the  ineident  wave  height 
representing  a  mixing  length  and  a  a  ealibration  faetor  of  0(1).  Negleeting  wave  refleetion  or  transfer 
to  lower  frequeney  waves,  all  incident  wave  energy  is  dissipated  by  wave  breaking  within  the 
surfzone,  and  the  surf  zone-averaged  wave  energy  dissipation,  ,  can  be  obtained  from  the  wave 
energy  balance  by  integrating  over  the  surfzone  width,  X5,  yielding: 


1  ? 

Xg  0  ^  s 


(2) 


Outside  the  surf  zone,  the  incident  energy  flux  is  constant  and  can  therefore  be  evaluated  with  the 
measurements  obtained  at  13  m  depth  (Table  1)  whereas  the  energy  flux  at  the  shore  line,  x  =  Xs , 
equals  zero.  Hence,  the  surfzone  averaged  wave  energy  dissipation  is  calculated  : 


K) 


S  ' x\h=-\3)m 


(3) 


with  the  incident  wave  energy  given  by: 


16 


(4) 


and  the  group  velocity  obtained  from  linear  wave  theory: 
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^71 

where  co  =  —  with  Tp  as  the  period  of  the  peak  of  the  wave  spectrum  and  the  wave  length  k  is 

'^p 

obtained  from  the  linear  dispersion  relation: 

a/  =  gktanhk/z 

(6) 


The  cluster-drifter  based  analysis  of  the  turbulent  eddy  diffusivity,  k^^ ,  represents  a  surf  zone-averaged 

quantity  (90%  of  the  drifters  stay  within  the  surfzone).  Using  the  values  for  waves  in  13  m  water 
depth  and  an  optimized  constant  value  of  a  =  0.7  and  X5  to  125  m,  there  is  a  significant  correlation 
between  the  estimated  surfzone  averaged  wave-breaking  induced  turbulent  eddy  viscosity,  ,  and  the 

estimated  turbulent  eddy  diffusivity  k^  at  the  95%  confidence  interval  (R  =  0.95). 
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Figure  3.  absolute  diffusivity  (one-particle)  statistics  (Kxx  (left)  (right)  plotted  with  sampling  error 
bars  S(t)  [dotted  lines]  against  time  for  yeardays  124  (top  row)  and  139  (bottom  row).  Peaks  in  Kxx 
occur  at  increments  of  ~300  s  corresponding  to  the  oscillations  in  the  auto-covariance  functions 
and  the  circulation  times  of  the  rip  circulation  cells  [MacMahan  et  al.  2008a].  The  quick  growth  of 
Kxx  for  F<100  s  followed  by  a  sharp  decrease  is  reflected  in  early  rapid  growth  in  Cxx  followed  by 
slower  expansion.  At  the  time  of  the  initial  peak  in  Kxx,  the  cross-shore  patch  size  Cxx  ~  25  m, 
indicating  that  for  a  basically  Gaussian  distribution,  96%  of  the  drifters  are  contained  within  4a xx  ~ 
100  m,  roughly  the  cross-shore  rip  current  dimension  (within  15%  difference),  indicating  the 
particles  diffuse  quickly  until  they  span  the  entire  surf  zone  at  which  time  cross-shore  diffusivity 
slows  as  particles  oscillate  between  the  seaward  edge  of  the  surf  zone  and  the  shoreline  [MacMahan 
et  al,  2008a].  On  yearday  139  the  cross-shore  absolute  diffusivity,  Kxx,  peaks  at  2.8  m^/s  then  drops 
to  near  zero  as  cross-shore  patch  expansion  stops  almost  completely  after  the  initial  dispersion.  The 
difference  in  diffusivity  is  associated  with  the  difference  in  the  flow  patterns. 


Delft3D  Model  Calculations 

A  3D  wave  and  flow  model  (DelftSD)  has  been  used  to  compute  the  wave  and  velocity  field  at  the 
wave-group  time  scale.  Comparison  with  in-situ  ADCP  measurements  show  that  the  hourly-mean 
variation  is  qualitatively  well  represented  (Figure  4).  Important  differences  do  occur,  most  notably  in 
the  prediction  of  the  vertical  structure  of  the  rip  current  flow  outside  the  outer  surf  zone  where  the 
observed  strong  vertical  attenuation  is  not  present  in  the  computations. 
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Figure  4  Hourly  mean  flow  velocity  at  ADCP4  (outer  surfzone).  Upper  panels  computed  cross-shore 
(left)  and  alongshore  (right)  Eulerian  velocity  (positive  offshore  and  downcoast).  Lower  panels: 

similar  for  observed  velocities. 
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Figure  5  Left  panel:  GPS-inferred  mean  velocities  for  the  deployment  on  yearday  124  compared 
with  computed  mean  surface  velocities  including  Stokes  drift  (GLM-velocities  indicated  in  yellow) 
and  without  (Eulerian  velocities  indicated  in  white).  Middle  panel:  Detailed  view  of  both  GLM  and 
Eulerian  mean  surface  flow  velocities  showing  the  stronger  curvature  for  the  GLM  velocities. 


Averaging  the  wave-group  resolving  surface  velocities  over  the  duration  of  the  drifter  deployment  on 
yearday  124  yields  the  mean  surface  velocity  field  which  is  compared  with  the  GPS-drifter  inferred 
velocities  (Figure  5).  The  computations  show  similar  flow  circulations  as  the  observations  with 
onshore  velocities  over  the  shallow  shoals  feeding  into  the  rip  currents.  Details  on  the  flow  divergence 
and  convergence  of  the  asymmetric  flow  circulations  are  also  well  represented.  The  offshore  extend  of 
the  rips  is  limited  due  to  the  fact  that  the  offshore  flow  quickly  diverges  resulting  in  mean  vortical 
circulations.  The  GLM  (including  Stokes  drift)  and  Eulerian  (without  Stokes  drift)  mean  flow 
velocities  are  very  similar  in  a  qualitative  sense  (Figure  5).  Flowever,  a  detailed  look  of  the  flow 
velocities  at  the  outer  surf  zone  shows  that  the  Eulerian  flow  field  generally  has  stronger  offshore 
directed  velocities  in  the  rip  channels  (right  panel  in  Figure  5).  The  alongshore  velocity  are  only 
moderately  affected  and  as  a  result  the  GEM  velocities  exhibit  a  stronger  rotation  (i.e.  smaller  radius) 
compared  with  the  Eulerian  velocities. 
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Figure  6  Left  panel:  Initial  drifter  locations  (black  circles)  and  calculated  trajectories  (yellow  dots) 
using  the  GLM  velocities  foryearday  124  for  hour  10.  Right  panel:  Similar  for  Eulerian  velocities. 
Position  of  the  outer  surfzone  edge  is  indicated  by  red  dashed  line.  Bathymetry  and  depth  contours 

in  m  given  as  a  reference. 


The  calculated  surface  velocity  fields,  at  the  wave  group  time  scale,  are  used  to  predict  the  trajectories 
of  virtual  surface  drifters.  This  is  done  for  GLM  velocities  (including  Stokes  drift)  and  Eulerian 
velocities  (without  Stokes  drift)  (Figure  6),  where  the  latter  result  generally  in  a  nearcomplete  exodus 
of  the  drifters  within  the  hour.  Using  the  GLM  velocities  the  majority  of  the  drifters  are  still  within  the 
surf  zone  after  one  hour.  This  is  a  result  of  the  fact  that  the  GLM  velocities  exhibit  a  stronger  rotation 
compared  with  the  Eulerian  velocities  (Figure  5).  The  GEM  derived  results  are  consistent  with  the 
field  observations,  whereas  the  Eulerian  results  over-predict  the  number  of  drifter  exits  by  an  order  of 
magnitude.  This  stresses  the  importance  of  including  Stokes  drift  in  calculating  the  transport  of  surface 
floating  material. 
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Figure  7  Estimated  percentage  of  the  mean  hourly  GLM  (solid  red  line)  and  mean  Eulerian  (solid 
blue  line)  drifter  exists  for  yeardays  124  through  126  as  function  of  the  tidal  elevation  plus  or  minus 
one  standard  deviation  (denoted  by  the  vertical  bars).  Observations  of  percentage  drifter  exits  as 
function  of  tidal  elevation  indicated  by  green  dots  on  yeardays  124, 125, 127  and  130  corresponding 

to  times  of  normally  incident  waves  during  RCEX. 


The  resulting  tidally-averaged  retention  rates,  eorresponding  to  the  pereentage  of  drifters  retained  in 
the  surf  zone  within  one  hour  are  0.82  for  the  GLM  velocities  and  0.15  for  the  Eulerian  velocities 
(Figure  7).  This  corresponds  to  retention  relaxation  time  scale  of  0(5)  hours  for  the  GLM  velocities 
and  less  than  1  hour  for  the  Eulerian  velocities.  As  a  result  the  retention  time  of  floating  surf  zone 
material  on  a  rip  channeled  beach  is  significantly  longer  than  originally  anticipated.  This  also  holds  for 
other  suspended  material  such  as  sediment  fines,  algae,  bubbles  and  other  organic  matter  tracking  the 
fluid  motions,and  are  therefore  subject  to  Stokes  drift.  However,  their  retention  also  depends  on  the 
vertical  advection  given  the  fact  these  agents  are  not  confined  to  the  water  surface  only  and  the  GLM 
velocities  vary  with  depth. 

In  summary,  the  model  calculations  closely  match  the  rip-current  surface  flow  field  inferred  from 
GPS-equipped  drifter  trajectories.  Also  the  ADCP-observed  temporal  velocity  variability  over  a  period 
of  three  days  with  normally  incident  swell  is  reproduced  by  the  computations,  although  the  strong 
vertical  attenuation  of  the  sub-surface  rip  current  velocities  at  the  most  offshore  location  is  not  well 
predicted.  Seeding  the  surf  zone  with  a  large  number  of  equally  spaced  virtual  drifters  the  computed 
instantaneous  surface  velocity  fields  are  used  to  calculate  the  drifter  trajectories  for  the  duration  of  an 
hour.  Surfzone  retention  is  examined  by  comparing  the  number  of  drifters  that  exit  the  surf  zone 
compared  to  the  total  number  of  active  drifters  during  one  hour.  Collecting  the  hourly  drifter  exits, 
good  agreement  with  the  observed  surfzone  retention  is  obtained  provided  that  Stokes  drift  is 
accounted  for  in  the  modeling  of  the  computed  drifter  trajectories.  By  including  Stokes  drift  O(20)\% 


12 


of  the  drifters  exit  the  surf  zone  within  one  hour.  Without  Stokes  drift,  the  estimated  number  of  drifter- 
exits  is  0(85)%,  an  order  of  magnitude  larger  than  the  0(10)%  of  exits  for  the  field  observations. 

IMPACT/APPLICATIONS 

Understanding  Lagrangian  rip  eurrent  behavior  is  eritieal  for  beaeh  safety,  water  quality,  water  elarity, 
nearshore  mixing,  and  sediment,  biologieal,  and  nutrient  exehange.  The  detailed  spatial  maps  of 
veloeity  and  vortiey  will  inerease  our  understanding  of  rip  eurrent  dynamies.  These  new  observations 
are  evalauted  with  DelftSD.  The  inelusion  of  Stokes  in  the  prediction  of  the  fate  of  floating  surface 
material  is  essential.  This  finding  has  important  consequences  for  calculating  the  transport  of  surface 
drifters  but  also  sediment  fines,  algae,  bubbles  and  other  organic  matter  and  as  a  result  will  also  impact 
the  predictions  of  both  water  clarity  and  water  qualitity. 
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